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ABSTRACT 


This research investigates thermal buckling characteristics of rectangular panels subjected to different types of 
humped temperature profile heating. Minimum potential energy and finite-element methods are used to calculate 
the panel buckling temperatures. The two methods give fairly close thermal buckling solutions. “Buckling 
temperature magnification factor of the first kind, r]” is established for the fixed panel edges to scale up the 
buckling solution of uniform temperature loading case to give the buckling solution of the humped temperature 
profile loading cases. Also, “buckling temperature magnification factor of the second kind, is established for the 
free panel edges to scale up the buckling solution of humped temperature profile loading cases with unheated 
boundary heat sinks to give the buckling solutions when the boundary heat sinks are heated up. 

NOMENCLATURE 

Acronyms 


BC 

JLOC 

SPAR 

TPS 


AR 

a 

V 

^mnkl 

Bmn 

b 

Cmn 

DiJ 


Dqj, Dqy 


boundary condition 
joint location 

Structural Performance and Resizing 
thermal protection system 


Symbols 

Fourier coefficient of trial function for out-of-plane displacement w, in. 


extensional stiffness of orthotropic plates, = - 


t.vE _ 

^21 - ~ ’ ^22 - 


1 - vv^. 


A^ 


66 = 


= a/b, plate aspect ratio 
length of plate, in. 

coefficient of characteristic equation 
Fourier coefficient of trial function for , in/in. 
width of plate, in. 

Fourier coefficient of trial function for y , in/in. 

El 

bending stiffness of plate, D, 1 = —-— >--ti - 

1 — vv 1 — vv \ — V , 

xy yx xy yx xy yx 

E I 

^66 = in-lb 


transverse shear stiffness in xz-, yz-planes, Ib/in. 



E 


E.,Ey 



a, a^, Oy 

a^y 

rxv yyz 

7] 

K 

K 


Young’s modulus, Ib/in^ 

Young’s modulus, Ib/in^ 
shear modulus, Ib/in^ 

moment of inertia, per unit width, with respect to plate centroidal axis, I^ ’ in^/in. 

index, 1, 2, 3,.... 
index, 1, 2, 3,.... 
index, 1, 2, 3,.... 
index, 1, 2, 3,.... 

index, 1, 2, 3,...., number of buckle half waves in x-direction 
thermal forces, Ib/in. 

index,l, 2, 3,...., number of buckle half waves in y-direction 
temperature (measured from room temperature), °F 
temperature for constant temperature profile heating, °F 
material temperature 

peak temperature of dome-shaped (or roof-shaped) temperature profile heating, °F 

boundary heat sink temperature, °F 

thickness of plate, in. 

x-displacement, in. 

y-displacement, in. 

out-of-plane displacement, in. 

rectangular Cartesian coordinates 

coefficient of thermal expansion, in/in-°F 

coefficient of thermal shear distortion, in/in-°F 

transverse shear strain in xz- and yz-planes, in/in. 

numerical factor in buckling equation written for specified panel edge condition 
= (, “buckling temperature magnification factor of the first kind” = buckling 
temperature dome temperature loading case (fixed or free support condition) 

divided by the buckling temperature of uniform temperature loading case (fixed 

boundaries) 

thermal buckling eigenvalue associated with constant temperature profile heating 
thermal buckling eigenvalue associated with dome-shaped (or roof-shaped) profile heating 
Poisson ratio 


2 



§ = g) _Q “buckling temperature magnification factor of the second 

kind” = buckling temperature of dome temperature loading case with any heat sink 

temperature (free support condition) divided by the buckling temperature {T of 

dome temperature loading case with zero heat sink temperature = 0 (free support 
condition) 

p density, Ib/in^ 

Oj tensile strength, Ib/in^ 

stress in x-direction, Ib/in^ 

Oy yield stress, Ib/in^ 

Oy stress in y-direction, Ib/in^ 

T shear stress, Ib/in^ 

( )g,. critical value at buckling 


INTRODUCTION 

Hypersonic aircraft are subjected to severe aerodynamic heating during flights. To maintain the structural 
integrity under high temperature environment, the vehicle structural design concepts of hypersonic flight vehicles 
are different from those of low Mach number aircraft. Depending on the operating temperature range, vehicle 
structures may be called “hot” structures or “warm” structures. The “hot” structures fabricated with high 
temperature alloys can operate at elevated temperatures in the range of 1000 °F to 1500 °F. If fabricated with the 
carbon/carbon composite material, the operating temperature of the “hot” structures could go as high as 3000 °F. 
Typical candidate “hot” structural components for hypersonic flight vehicles are hat-stiffened panels fabricated 
with either monolithic titanium alloys or metal-matrix composite materials; honeycomb sandwich panels 
fabricated with super-alloy; and carbon/carbon composite eleven (or body flap). 

The “warm” structures are fabricated with lightweight materials, such as aluminum, and can operate only up to 
moderate temperature limit of 350 °F. The space shuttle orbiter is a good example of the “warm” structure. The 
entire vehicle is protected with the thermal protection system (TPS) to shield the aluminum substructure from 
overheating beyond the “warm” temperature limit. 

Hot structural panels are usually fastened to the cooler substructures that function as heat sinks because of less 
heating. Thus, even under the uniform surface heat flux, the temperature distribution over the hot structural panels 
will not be uniform but looks like camel-humped shape (refs. 1-3). This camel-humped shaped temperature 
distribution is always observed in supported hot structural panels. The buckling behavior of the panel with heat 
sinks is therefore quite different from the case without heat sinks. Earlier Ko extensively studied thermal buckling 
problems of hot structural panels under uniform temperature profile heating (without heat sinks effects) 
(refs. 4-10). 

This report studies the thermal buckling behavior of rectangular panels under different types of 
humped-shaped temperature profile heating to simulate the existence of the supporting boundary heat sink and also 
studies the effect of the heat sink temperature on the panel buckling temperature. The results are compared to the 
results of uniform temperature profile heating cases without the heat sink effect. 
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HUMPED TEMPERATURE PROFILES 


Under uniform surface heating, because of the existence of the supporting cooler boundary heat sinks, 
experimental data have shown that the temperature distribution over a hot structural panel is usually 
flattened-dome shaped (frg. 1, ref. 1). For performing thermostructural analysis, the actual temperature profile may 
be idealized with simple mathematical surfaces such as dome-shaped or roof-shaped temperature profiles 
described in the following sections. 

Dome-Shaped Temperature Profile 

The dome-shaped temperature profile for a rectangular panel (length a, width b, thickness t^) is described by 
the following mathematical function: 


T{x,y) = r +(r -rjsin —sin^ 
^ ^ a b 


( 1 ) 


where is the peak temperature at the panel center and is the boundary heat sink temperature (fig. 2). 
Equation (1) gives a “sine-sine” surface elevated upward by an amount (fig. 2). Keep in mind that the 
temperatures in this report are the temperature differentials above the room temperature. 

Equation (1) is the first-order approximation of the actual temperature profile (fig. 1). Eor more accurate 
representation of the actual temperature profile, double Eourier series representation is necessary; however, this 
representation is not the scope of the present report. 

When the boundary temperature goes to zero (fig. 3), equation (1) becomes: 

T{x,y) = r sin —sin^ (2) 

° a b 


In order to obtain closed-form buckling solutions by means of the minimum potential energy theory, the functional 
form given by equation (2) is used for dome-shaped temperature loading. 

Eor a special case of uniform temperature profile heating (fig. 4), equation (1) degenerates into: 

T{x,y) = T^ = T^ = (3) 

Thermal buckling problems of rectangular sandwich panels under the constant temperature profile heating 
[eq. (3), fig. 4] was analyzed earlier by Ko (refs. 4—7) using the minimum potential energy theory and 
finite-element method. 


Roof-Shaped Temperature Profile 

The actual dome-shaped temperature profile (fig. 1) may also be approximated by a simple temperature 
distribution surface called roof-shaped temperature profile (peak temperature at the profile central flat zone and 
boundary heat sink temperature T^) (figs. 5 and 6). When = T^, the roof-shaped temperature profile also 
degenerates into uniform temperature profile (fig. 4). 
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PROBLEMS INVESTIGATED 


This report investigated thermal huekling behavior of reetangular panels (length a, width b, thiekness t^) 
heated under the dome-shaped (figs. 2 and 3) and roof-shaped (figs. 5 and 6) temperature profiles. The main 
aspeets of this researeh are: 

1. To study the effeet of the dome-shaped (or roof-shaped) temperature profile on the huekling behavior of the 
reetangular panels with fixed edges. Then, the study will establish the “huekling temperature magnifieation 
faetor of the first kind” to seale up the huekling solution of uniform temperature profile heating ease to give 
the huekling solution of the dome-/roof-shaped temperature profile heating ease. 

2. To study the effeet of the heat sink temperatures on the huekling temperatures of reetangular panels with 
free edges under dome-shaped (or roof-shaped) temperature profile heating. Then, the study will establish 
the “huekling temperature magnifieation faetor of the seeond kind” to seale up the huekling solution of 
dome (or roof) temperature loading ease with unheated boundary heat sinks to give the huekling solutions 
when the boundary heat sinks are heated up. 


FINITE-ELEMENT BUCKLING ANALYSIS 


This researeh will use the Struetural Performanee and Resizing (SPAR) finite-element eomputer program (ref 
11) to eonduet the linear elastie thermal huekling analysis of the reetangular panels. 

Panel Geometry 

All the reetangular panels analyzed have thiekness = 0.09 in. with the different dimensions in length a and 
width b listed in Table 1. 


Table 1. 

Dimensions of reetangular panels; 

= 0.09 in. 

a, ir 

1 . b, in. 

alb 

6 

6 

1.0 

9 

6 

1.5 

12 

6 

2.0 


Figure 7 shows the quarter model generated for the reetangular panel. The quarter model has 625 joint 
loeations (JLOCs) and 576 E43 four-node plate elements. 

Material Properties 

The reetangular panels are made of high temperature alloy Haynes® 230® (Haynes International, Koyomo, 
Indiana) with the room temperature properties listed in Table 2. 
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Table 2. Material properties of Haynes 230 alloy at room temperature. 


E 

30.6 X 10^ Ib/in^ 

V 

0.31 

P 

0.324 Ib/in^ 

a 

7.0 X 10"^ in/in-°F 


Haynes 230 high temperature alloy was used in the fabrication of Hyper-X (X-43) hypersonic research vehicle 
for Mach 7-10 mission (ref. 3). Appendix A gives the temperature-dependent material properties of Haynes 230 
alloy. 


Boundary Conditions 

The following different panel boundary conditions were considered in the finite-element analysis: 

1. 4S fixed — four edges simply supported, with no in-plane displacements [fig. 8(a)]. 

a. u{x,y) ^ 0,v{x,y) 0 everywhere except boundaries. 

b. u{x,y) = 0, v{x,y) = 0 everywhere in the panel. 

2. 4S free — four edges simply supported and slides freely along the lubricated guides, which can have free 
in-plane motions [fig. 8(b)]. 

3. 4C fixed — four edges clamped, with no in-plane displacements [fig. 9(a)]. 

a. u{x,y) 7^ 0, v(x,y) 0 everywhere except boundaries. 

b. u{x,y) = 0, v(x,y) = 0 everywhere in the panel. 

4. 4C free — four edges clamped and slides freely along the lubricated clamping guides, which can have free 
in-plane motions [fig. 9(b)]. 


Thermal Loads 

The thermal load inputs used for the dome (or roof) temperature loading and constant temperature loading are 
as follows. 

1. Dome-Shaped (or Roof-Shaped) Temperature Profile Heating 

For the dome-shaped [eqs. (1) and (2), figs. 2 and 3] and roof-shaped (figs. 5 and 6) temperature load inputs to 
the SPAR program for the eigenvalue calculations, the unit peak temperature = 1 °F was used. The heat sink 
temperature was allowed to vary over the range T^/T^= 0~1 (i.e., T^/T^ = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.8, 0.9, 1.0). Then the eigenvalue (scaling factor) calculated from SPAR program will give the buckling 
temperature namely: 


(r ) =Axl=A^ 

^ o'er o o 


(4) 
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If the input peak temperature Tg is doubled (i.e., = 2 °F), then the ealeulated eigenvalue will be redueed to half 

of Xg given in equation (4) beeause the produet {TgX Xg) remains eonstant. 

The reason for using the whole range of T/Tg = 0~1 is to study the effeet of the heat sink temperature on 
the panel buekling temperature {Tg)^^ when the panel edges are free to move. Keep in mind that the panel 
buekling temperature {Tg)^^ will inerease with the inereasing heat sink temperature beeause of relaxation of 
thermal expansion eonstraint exerted on the panel by the boundary heat sink. 

2. Uniform Temperature Profile Heating 

As mentioned earlier, the uniform-temperature loading ease (fig. 4) is a degenerative ease of dome-shaped (or 
roof-shaped) temperature profile heating. Temperature load of = 1 °F is ehosen as input to all nodes of the 
finite-element model so that the eigenvalue X^ ealeulated from SPAR program will give the buekling temperature 
namely: 


This buekling solution is for the panels with fixed boundary supports. When the panel boundaries ean have free 
in-plane motions (free supports), the eonstant temperature loading ease obviously ean never induee thermal 
buekling (i.e., buekling temperature goes to infinity). 

3. Material Property Iterations 

In the ealeulations of buekling temperatures, material property iterations are required to obtain aeeurate 
buekling temperatures. Namely, the material properties at eertain assumed material temperature must be 
updated until matehes the ealeulated buekling temperature {Tc\r' 12 diseusses in detail this 

iteration proeess. Usually it takes only two to three material iterations to yield aeeurate buekling temperatures. 

BUCKLING TEMPERATURE MAGNIFICATION FACTORS 


As diseussed later in this report, for the fixed supports (4S or 4C) [figs. 8(a) or 9(a)], the buekling temperature 
{To)cr ll*® dome (or roof) temperature loading ease is mueh higher than the buekling temperature {Tg)g^ of 
the uniform temperature loading ease. Therefore, a “buekling temperature magnifieation faetor of the first kind, r]” 
defined as: 


To)cr 


will be used to indieate how many times the buekling temperature {Tg)^^ of the dome (or roof) temperature 
loading ease (fixed supports) is magnified from the buekling temperature (T^)^^ of the uniform temperature 
loading ease with fixed supports (fundamental ease). 

In reality, the hot panel attempts to expand under heating; but its expansion is resisted by the eooler boundary 
substruetures (heat sinks) that expand less. This boundary eonstraint is the eause of thermal buekling of the panels. 
Sueh substrueture eonstraints will gradually relax as the substruetures are heated up, resulting in higher panel 
buekling temperature. In order to diseover how the heat sink temperature affeets the panel buekling temperature 
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{To)cr fo*" dome (or roof) temperature loading ease with free supports, another “huekling temperature 
magnifieation faetor of the seeond kind, defined as: 



will he used to indieate how many times the huekling temperature for the dome (or roof) temperature 

loading ease with free support (4S or 4C) is magnified when the heat sink temperature inereases from = 0 
(no heat sink thermal expansion) to eertain nonzero value T^ 7 ^ 0 (with heat sink thermal expansion). 


MINIMUM POTENTIAL ENERGY BUCKLING ANALYSIS 


When the [u, v} displaeements are zero everywhere in the plate ineluding the boundaries (4S fixed-h and 4C 
fixed-h eases), the dome-shaped temperature profile heating [eq. (2), fig. 3] will generate identieal dome-shaped 
biaxial stress distributions with no shear stresses. Beeause of this ideal eondition, it is possible to obtain 
elosed-form huekling solutions for this partieular ease using the minimum potential energy theory. 

Energy Equations 

Strain energy assoeiated with bending Vj may be written as (refs. 5 and 6): 






( 8 ) 




Lo*y\o*y 2 Idyxdx dx\dy 

and the strain energy assoeiated with thermal loading V 2 is given by: 

1 / 


The dome-shaped temperature loading funetions may be expressed as follows (refs. 5 and 6): 


( 9 ) 


^^{x,y) = {A^^a^ + Ai2CCy)T — ( 10 ) 

Ny{x,y) = {A2ia^ + A22ay)T^sm^sm^ ( 11 ) 

Nly{x,y) = 0 ( 12 ) 

The preeeding thermal loading funetions (10)-(12) hold only for the 4S fixed-b and 4C fixed-b eases 
[i.e., M(x,y) = 0, v(x,y) = 0 everywhere]. 




For other fixed panel eases under dome temperature heating, the in^lane displaeements of panel interior points 
are nonzero [i.e., u{x,y) 0, v{x,y) 0], and the thermal loads , N} are eomplieated funetions of [x,y}; 

and thus it will he too eumhersome to apply the minimum potential energy method. Therefore, SPAR was used to 
obtain quiek huekling solutions. 


Deformation Functions 

The deformation funetions for the 4S fixed-h ease and 4C fixed-h ease may he expressed in the following 
forms (refs. 5 and 6). 

Case 1. Four edges simply supported (4S fixed-h ease) 


w(x,y) =22 

y^^(x,y) =22 

m = In = 1 

= 2 2 


n = In = 1 

Case 2. Four edges elamped (4C fixed-h ease) 


( 13 ) 


(14) 


( 15 ) 


w(x,y) = sin^sin^ V V Asin^sin^ 

a bZjZj a 


(16) 


y, (x,y) = eos^sin^ V V R sin^^l^sin^^ 
a b Zj Zj ” a b 


+ sin^sin^ y y ruBcos'^sin^ 

a b Lj Lj a b 


( 17 ) 


/ N ■ Jix ■ mjix ■ njty 

,,Jx,y) = sm — eos^ \ \ C sm-sm—^ 

a b Zj Zj a b 

■ Jtx ■ Jiy ^ ■ niJix njiy 

+ sm — sm^ \ \ nC sm-eos —^ 

a bZjZj a b 


(18) 


Buckling Equations 

Applieation of the minimum potential energy theory yields the following homogeneous simultaneous 
eharaeteristie equations (refs. 5 and 6) for eaeh set of integral values {m,n} (or mode shape) written for orthotropie 
reetangular panels under dome-shaped temperature profile heating [eq. (2)]: 
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(19) 


where the stiffness/geometry parameter is given hy: 


^mnkl 


1 11 
^ ‘^mnkl 


^mnkA^mnkl^mnkl ~ ^mnkfimnkO ^mnkA^mnkl^mnkl ~ ^mnkl^mnW 
22 33 23 32 

^mnkl^mnkl ~ ^mnkl^mnkl 


( 20 ) 


and where are the eoeffieients of eharaeteristie equations. 

Appendix B gives funetional expressions of whieh is taken from Appendix C of referenee 6 (or 

Appendix A of ref. 7) with the loading terms eontaining {k^, ky] removed. 

For the simply-supported edges (4S fixed-h ease), the value of ^ in equation (20) is ^ = 32; and the thermal 
foreing term P^nnkl equation (19) has the following form (speeial integrals given in Appendix C were used in the 
ealeulations). 


n l/T T Arn^ + k^)-l ^ (n^-i-/^)-l 

Pmnkl = - (^ll«;c + '^12«P-3-+ (^21 «;c +'^22«P-^- 


mnkl 


{{m + k) -l][(m-k) -l][(n + 0 -!][(«-/) -1] 


( 21 ) 


m ± k = even, n ± I = even 

For the elamped edges (4C fixed-h ease), the value of ^ in equation (20) is ^ = 24, and P^nkl equation (19) 
takes on the following form (speeial integrals given in Appendix C were used in the ealeulations). 


n mnkL . . , 

Pmnkl = ^(^ll«x + ^12^) 


Am + k)^ - l][(m - kY - 1] \_{m + k)^ - 9][(m - k)^ - 9] 

_ =1 _+_^_ 

{{n + lf-\ ][{n -if-l] [{n + if -9][{n - if - 9] 


.2 


.2 


.2 


( 22 ) 


2 ('^2l“x''''^22“>’- 


[{n + l) -l][(n-/) -1] [{n + l) -9][{n-l) -9] 
-1 , 1 


[{m + kf-l ][{m -kf-l] [{m + kf - 9][{m - kf - 9] J 

m ± k = even, n ±l = even 
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For the unknown deflection coefficients in equation (19) to have nontrivial solution, the determinant of the 
coefficients of unknown Aj^i of simultaneous equations written out from equation (19) must he set to zero for 
eigenvalue extractions. The largest eigenvalue (l/T^) will then give the lowest critical buckling temperature 

Appendix D shows the determinants of the coefficients of simultaneous equation (called buckling equations) 
written out from equation (19) up to order 12 for the cases m±n = even (symmetrical buckling) and m±n = odd 
(antisymmetrical buckling) for the 4S and 4C cases. The symmetrical buckling case (m ± n = even) gives the 
lowest buckling temperature associated with the lowest buckling mode. 

Buckling Temperatures for Isotropic Cases 

The purpose of the present buckling analysis is to study the effect of dome-shaped temperature profile heating 
on the buckling of the rectangular panels and not the effect of material anisotropy. Therefore, only the isotropic 
panels are considered. 

1. Simply-Supported Edges 

For the simply-supported case (4S fixed-b case) ignoring the transverse shear effect, [eq. (20)] and Pm^kl 
[eq. (21)] take on the following forms: 


Et"jt 


384(1 -v )l a E 


(23) 


Et^a (m +k)-l {n -l-/)-l 
1 - V ^2 2 


(24) 


mnkl 


[{m + kt-\ ][(m -kt-\ ][{n + lt-\ ][(u -/)-!] 


m±k = even, n±l = even 

The first-order solution (m = u = k = / = 1) of the buckling temperature for the dome-shaped heating may 
be calculated through combining equations (19), (23), and (24) by setting m = n = k=l=l. 




3;r^^ [l l' 

128«(1 + v )[ 72-^-2 


(25) 


For the constant temperature heating case, the first-order (m = \, n = V) buckling temperature (7’(,)cr for the 
isotropic panel can be calculated from equations (28)-(30) of reference 7. With the transverse shear effect 
neglected, the first-order buckling temperature (7’c)cr for this case has the following form: 
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(26) 




12a(l + v)l 2 


Dividing equation (25) by equation (26) yields the buekling temperature magnifieation faetor of the first kind 
Tj for the simply-supported first-order ease (m = u = k = / = 1) in a neat form: 




9;r 

32 


(27) 


that is independent of panel aspeet ratio. 

For the 4S ease, the values of rj based on the seeond and third-order buekling solutions (whieh depend on the 
aspeet ratio alb) were also obtained. The RESULTS seetion presents these values. 

2. Clamped Edges 

For the elamped ease (4C fixed-b ease), the first-order solution (m = u = k = Z = 1) for the “dome-shaped” 
temperature ease is: 


(28) 


675(1 -v)l 2 ^2 


(29) 


The first-order solution (m = u = k = / = 1) of the buekling temperature for the elamped panel under the 
dome-shaped heating may be ealeulated by eombining equations (19), (28), and (29) to yield: 




64^a(l-i-v) 


(30) 


The first-order solution (m = n = k = I = 1) for the eonstant temperature ease ealeulated from equations 
(28)-(30) of referenee 7 gives: 




3 2 3 

— -I---I- — I 

.4 2,2 ,4 

Jt t \a a b b 


9a(l -I- v) 


2 ,2 


( 31 ) 
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The ratio of equations (30) and (31) gives the huekling temperature magnifieation faetor of the first kind f] for 
the 4C ease first-order solution (m = n = k = I = 1) in a neat elosed form as: 


r] = 




= 1.6265 
64 ^ 


(32) 


whieh, similar to the 4S ease [eq. (27)], is also independent of panel aspeet ratio. As diseussed in a following 
seetion, for the 4C ease, the first-order huekling solution gives a quite aeeurate value of rj. Therefore, higher order 
values of r] were not ealeulated. 


RESULTS 

The following seetions present the results of the thermal huekling analyses of the reetangular plates under 
different heating profiles. 


Buckling Temperature Magnification Factors 

Figures 10 and 11 show the huekled shapes of the isotropie square panel {alb =1) under different temperature 
profile loading and under different boundary eonditions. Notiee that the panel huekling shapes are insensitive to 
the temperature loading funetions for eaeh set of edge eonditions. 

Table 3 lists the huekling temperature magnifieation faetor of the first kind 17 = ealeulated 

from SPAR program for the isotropie reetangular panels with different aspeet ratios under dome-shaped 
temperature profile heating with zero heat sink temperature [T^ = 0 , eq. ( 2 )]. 


Table 3. Buekling temperature magnifieation faetor r] 
for reetangular plates under dome-shaped temperature 
profile heating; = 0 ; fixed boundaries. 




rj 


BC \ alb 

1 

1.5 

2 

4S 

1.9721 

(2.0034) 

1.9603 

1.9453 

4C 

1.8035 

(1.8299) 

1.7848 

1.7320 


() Sandwich panel (ref 7). 


Notiee that r] is insensitive to the ehange of panel aspeet ratio alb and that the 4S eases have slightly higher 
values of rj than the 4C eases. 

In Table 3, the rj values for the square sandwieh panel (alb = 1) are also shown in the parentheses. Those 
values were ealeulated from the previous SPAR program developed for the buekling of sandwieh panel (ref. 7). 
Notiee the proximity of the f] values for solid and sandwieh panels. 
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Table 4 compares the values of rj calculated from the minimum potential energy theory and from the SPAR 
program for the special case of square panel {alb = 1) under the dome-shaped temperature profile heating = 0, 
eq. (2), frg. 3], and with constraint u = v = 0 everywhere. 


Table 4. Buckling temperature magnification factor rj for square plate {alb =1) 
under dome-shaped temperature profile heating; = 0\ u = v = 0 everywhere. 


BC Method 


7) 


Ko theory 

SPAR 

Percent difference 

4S 

2.4503*** 

2.5185 

2.71 


(2.6394)** 


4.80 


(2.7758)* 


10.22 

4C 

1.6265* 

1.5993 

1.70 


First-order solution. 
Second-order solution. 


*** Third-order solution. 

For the 4S case, the r] values were hand calculated up to the third-order buckling solutions. The third-order 
solution compares quite well with the SPAR solution with only 2.71 percent difference. Thus, higher order 
solutions were not pursued. 

For the 4C case, the first-order solution is extremely good with only 1.70 percent difference compared with the 
SPAR solution. Therefore, higher order solutions for f] were not calculated. 

Table 5 lists the buckling temperatures {Tg)^,. and the values of buckling temperature magnification factors 
[rj, calculated for the simply-supported (4S fixed/free) square panel {alb =1) under shifted dome-shaped 
temperature loading [eq. (1), fig. 2] under different heat sink temperatures 
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Table 5. Buckling temperatures of square plate {a/b = 1) under 
dome-shaped temperature profile heating; 4S fixed and free cases. 


TjTg 

4S Fixed 

4S Free 


(To)cn °F 

T] 

(Tglcr, °F 

§ 

0.0 

80 

2.00 

515 

1.00 

0.1 

73 

1.83 

568 

1.10 

0.2 

67 

1.68 

635 

1.23 

0.3 

62 

1.55 

718 

1.39 

0.4 

57 

1.43 

821 

1.59 

0.5 

54 

1.35 

960 

1.86 

0.6 

50 

1.25 

1,175 

2.28 

0.7 

A1 

1.18 

1,500 

2.91 

0.8 

45 

1.13 

2,070 

4.02 

0.9 

42 

1.05 

(4,213) 

8.18 

1.0* 

40 

1.00 

00 

00 


Uniform temperature loading. 

() Exceeded melting range 2,375 °F~2,500 °F of Haynes 230 alloy. 


Note that for the 4S fixed edges, the value of f] decreases with the increasing heat sink temperature T^. For the 
4S free edge case, ^ increases with the increasing heat sink temperature Ty When TJT^ = 1, the case degenerated 
into the uniform temperature loading case for which the buckling temperature goes to infinity for the 4S free edge 
case (i.e., thermal buckling does not occur). 

It must be mentioned that for the actual hot sandwich panel (ref. 1) TJT^ = 0.45 and for the HYPER-X wing 
panel (ref. 3) TJT^ = 0.54. 

Table 6 lists the buckling temperatures and the values of buckling temperature magnification factors 

[f], calculated for the clamped (4C fixed/free) square panel {alb =1) under shifted dome-shaped temperature 
loading [eq. (1), fig. 2] and under different heat sink temperatures Ty 
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Table 6. Buckling temperatures of square plate {alb = 1) under 
dome-shaped temperature profile heating; 4C fixed and free cases. 


TJTo 

4C Fixed 

4C Free 

(To)cn °F 

7] 

iTo)cr^ °F 

§ 

0.0 

194 

1.81 

810 

1.00 

0.1 

180 

1.68 

886 

1.09 

0.2 

167 

1.56 

982 

1.20 

0.3 

156 

1.46 

1,110 

1.34 

0.4 

147 

1.37 

1,272 

1.57 

0.5 

138 

1.29 

1,480 

1.83 

0.6 

130 

1.21 

1,770 

2.19 

0.7 

124 

1.16 

2,302 

2.84 

0.8 

118 

1.10 

(3,581) 

4.42 

0.9 

112 

1.05 

(6,907) 

8.53 

1.0* 

107 

1.00 

00 

00 


Uniform temperature profile heating. 

() Exceeded melting range 2,375 °F~2,500 °F of Haynes 230 alloy. 


Note that for the 4C fixed case, similar to the 4S fixed case, value of f] decreases with the increasing heat sink 
temperature T^. For the 4C free case, the values of ^ are very close to the corresponding ^ value of 4S free case 
(Table 5). Again, when reaches (i.e., TJT^ - 1), the case degenerated into the uniform temperature loading 
case for which the buckling temperature goes to infinity for the 4C free case (i.e., no thermal buckling). 

Table 7 lists the buckling temperatures and the values of buckling temperature magnification factors 

[r], calculated for the simply-supported (4S fixed/free) square panel {alb = 1) under shifted roof-shaped 
temperature loading (fig. 5) under different heat sink temperatures T^. 
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Table 7. Buckling temperatures of square plate {alb =1) under 
roof-shaped temperature profile heating; 4S fixed and free cases. 


TJTo 

4S Fixed 


4S Free 

iTo)cr^ °F 

7] 

iTo)cn °F 

§ 

0.0 

51 

.28 

622 

1.00 

0.1 

50 

.25 

690 

1.11 

0.2 

48 

.20 

767 

1.23 

0.3 

A1 

.18 

860 

1.38 

0.4 

46 

1.15 

985 

1.58 

0.5 

45 

1.13 

1,160 

1.86 

0.6 

44 

I.IO 

1,410 

2.27 

0.7 

43 

1.08 

1,764 

2.84 

0.8 

42 

1.05 

(2,587) 

4.16 

0.9 

41 

1.03 

(5,174) 

8.32 

1.0* 

40 

l.OO 

00 

00 


Uniform temperature profile heating. 

() Exceeded melting range 2,375 °F~2,500 °F of Haynes 230 alloy. 


Table 8 lists the buckling temperatures and the values of buckling temperature magnification factors 

[f], calculated for the clamped (4C fixed/free) square panel {alb = 1) under shifted roof-shaped temperature 
loading (fig. 5) under different heat sink temperatures T^. 


Table 8. Buckling temperatures of square plate {alb = 1) under 
roof-shaped temperature profile heating; 4C fixed and free cases. 


TJTo 

4C Fixed 

4C Free 

(Jo)cr^ °F 

7] 

(Jo)cr^ °F 


0.0 

132 

1.23 

1,125 

1.00 

0.1 

129 

1.21 

1,235 

1.10 

0.2 

126 

1.18 

1,370 

1.22 

0.3 

123 

1.15 

1,520 

1.35 

0.4 

121 

1.13 

1,720 

1.53 

0.5 

118 

1.10 

1,980 

1.76 

0.6 

116 

1.08 

(2,501) 

2.22 

0.7 

113 

1.06 

(3,335) 

2.96 

0.8 

111 

1.04 

(5,002) 

4.45 

0.9 

109 

1.02 

(10,004) 

8.89 

1.0* 

107 

1.00 

00 

00 


Uniform temperature loading. 


() Exceeded melting range 2,375 °F~2,500 °F of Haynes 230 alloy. 
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Figures 12 and 13, respectively, show the buckling temperatures (7’^)^^ (data from tables 5-8) plotted as 
functions of normalized heat sink temperature TJT^ for the fixed and free edge cases to illustrate entirely different 
buckling behavior of the fixed and free edge cases. Note from figure 12 that for the fixed edges, the buckling 
temperatures of the roof-shaped temperature profile case are much lower than the dome-shaped temperature profile 
case for each panel support condition. However, for the free edges (fig. 13), the reverse is true. 

Figure 14 shows the buckling temperature magnification factors [f], (data from tables 5-8) plotted as 
functions of normalized heat sink temperature TJT^ for easy visualization of the shapes of { r), curves. Note that 
Tj values for the roof case are much lower than the dome case for the same support condition. For the free edge 
cases, the ^ curves for different loading and edge support conditions stay very close. 

In-Plane Deformations 

When the panel is constrained only at the four boundaries, the panel interior points can have in-plane free 
motions if the temperature loading is nonuniform. Figure 15 shows the in-plane deformed shapes of the square 
panel under dome-shaped [fig. 15(a)] and roof-shaped [fig. 15(b)] temperature loading cases. As expected, the 
panel center regions expand more than the panel boundary regions as the result of heating temperature profiles. 
Figures 16 and 17, respectively, show the distributions of in-plane displacements u(x,y) and v(x,y) for the 
dome-shaped temperature loading. Note that {u, v] reach a maximum at points between the panel center and the 
boundaries. For the roof-shaped temperature profile case (figs. 18 and 19), the maximum points of [u, v] migrated 
toward the panel center. 


Thermal Stresses 

Figures 20 and 21, respectively, show the compressive thermal stresses (a^, o^} induced in a fixed square 
panel {alb = 1) under the unit dome-shaped temperature profile heating {T^ = 1 °F, = 0; eq. (2), fig. 3]. Even 

though the loading temperature profile is dome shaped, the distributions of thermal stresses [o^, o } are 
airplane-hanger shaped because u{x,y) 0 and v{x,y) 0 within the fixed boundaries. Figure 22 shows the 
distribution of shear stress which is wavy shaped. The peak magnitudes of the shear stress occur at the 
boundaries and near the panel comers. 

For the roof-shaped temperature profile heating (T^ = 1 °F, = 0) (fig. 6), the distributions of [ } are 

roughly trapezoidal hanger-shaped (figs. 23 and 24). The shear stress distribution (fig. 25) exhibits zero shear 
stress near the panel center region and a maximum shear near each panel comer. 

If the in-plane displacements [u, v] are constrained everywhere in the panel including the edges 
[i.e., u{x,y) = v(x,y) = 0], the resulting distributions of the biaxial stresses [ a^, o^} will be dome-shaped just like 
the input dome temperature distribution (figs. 26 and 27), and the in-plane shear stress diminishes {x^^ = 0). For 
the case of roof-shaped temperature profile (figs. 28 and 29), the distributions of the biaxial stresses [ , o^} 

reflect the input temperature profile. 

Table 9 lists the fundamental thermal stresses induced in the fixed square panel {a/b = 1) by the dome-shaped 
(or roof-shaped) temperature loading (7), = 1 °F) with zero heat sink temperature (7), = 0). 
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Table 9. Peak thermal stresses {CT^, Oy induced in square plate {alb = 1) under 
dome-shaped and roof-shaped temperature loads; Tg=l °F; = 0; fixed boundaries. 


Constraints of {w, v} 

Heating 

O^, Ib/in^ 

Oy Ib/in^ 

Ty.y Ib/in^ 

M = 0, V = 0, edges only 

Dome 

Roof 

-189.26 

-256.38 

-189.26 

-256.38 

31.27 

48.66 

M = 0, V = 0, everywhere 

Dome 

Roof 

-301.05 

-306.77 

-301.05 

-306.77 

0.00 

0.00 


Note that by setting u{x,y) = v(x,y) = 0 everywhere in the panel including the boundaries, the peak magnitudes 
of the biaxial stresses {o^, Oy] were raised respectively by 59 percent and 20 percent for the dome- and 
roof-shaped temperature profile cases. 


CONCLUDING REMARKS 

Thermal buckling characteristics of rectangular panels subjected to dome-shaped and roof-shaped temperature 
profile heating were investigated using the finite-element method and the minimum potential energy method. The 
key results include: 

1. “Buckling temperature magnification factor of the first kind, r]” was established to scale up the buckling 
solution of uniform temperature loading case to give the buckling solution of the dome-shaped (or 
roof-shaped) temperature loading cases. Also, “buckling temperature magnification factor of the second 
kind, was established to scale up the buckling solution of dome temperature loading case with unheated 
boundary heat sinks to give the buckling solutions when the boundary heat sinks are heated up. 

2. For the fixed boundary cases when the panel interior in-plane motions are not constrained (u^O, v 0 
except boundaries), the panel buckling temperatures under dome-shaped temperature loading are 
practically twice the buckling temperatures of the uniform temperature loading cases. 

3. For simply-supported case (4S fixed, u = v = 0 everywhere) under the dome-shaped temperature profile 
heating, the third-order solution of the “buckling temperature magnification factor of the first kind, r]” 
calculated from the minimum potential energy theory, agrees fairly well with the r] value calculated from 
the finite-element method. The solution difference is only 2.71 percent. 

4. For the clamped case (4C fixed, u = v = 0 everywhere) under dome-shaped temperature profile heating, the 
first-order solution of minimum potential energy gives quite accurate value of rj with onlyl.70 percent 
difference from the finite-element solution. 

5. For the fixed boundary case when the panel interior points can have free in-plane motions (w 0, v 0 
except boundaries), the distribution profiles of the compressive stresses [a^, Oy] under the dome-shaped 
(or roof-shaped) temperature loading have airplane-hanger-like shapes. 

6. For the fixed boundary case when the panel in-plane motions are constrained everywhere (w = v = 0 
everywhere), the distribution profiles of the compressive stresses {ct^, Oy] take on the dome shape (or roof 
shape) similar to the input dome-shaped (or roof-shaped) temperature profile. 

7. The “buckling temperature magnification factor of the first kind, r]” associated with the dome-shaped 
temperature profile heating is much greater than that for the roof-shaped temperature profile heating case 
for both 4S and 4C fixed edge conditions. 
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8. The “buckling temperature magnification factor of the second kind, associated with the dome-shaped 
and roof-shaped temperature profile heating are very close for 4S and 4C free edge conditions. 


Dry den Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, April 3, 2002 
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APPENDIX A 

TEMPERATURE-DEPENDENT MATERIAL PROPERTIES OF 
HAYNES 230 ALLOY 


Temperature-dependent material properties of Haynes 230; p = 0.324 Ib/in^. 


r, °F 

10^,lb/in2 

ay X 10^, lb/in2 

Oy X 10^, lb/in2 

a X 10 ^, in/in- 

°F V* 

70 

30.6 

125.4 

57.4 

7.0 

0.310 

200 

30.1 

(122.2) 

(55.0) 

7.1 

0.311 

400 

29.3 

(117.3) 

(51.3) 

12 

0.315 

600 

28.3 

(112.3) 

(47.7) 

lA 

0.318 

800 

27.3 

(107.4) 

(44.0) 

7.6 

0.321 

1000 

26.4 

102.5 

40.3 

7.9 

0.324 

1200 

25.3 

97.7 

39.5 

8.1 

0.330 

1400 

24.1 

87.7 

42.5 

8.3 

0.332 

1600 

23.1 

63.1 

37.3 

8.6 

0.334 

1800 

21.9 

35.2 

21.1 

8.9 

0.340 

2000 

20.7* 

19.5 

10.8 

9.2* 

0.343 


Estimated. 


() Interpolated. 
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APPENDIX B 

COEFFICIENTS OF CHARACTERISTIC EQUATIONS 

The characteristic coefficients appearing in equation (20) are defined in the following for different 

indicial and edge conditions (refs. 6, 7). 

Case 1: 4S condition 




(y) -'’-(y) 




^mnmn ^mnmn 


Case 2: 4C condition 


1. m = n=k = l = l 


>0 

= -[i20„Q’+4(D,; 

= -[l2022(f)’+4(D,; 
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■■it) ^Mt) 


2. m = k = l,n= 1^1 




Mt) 0 -^)] 


13 31 3 




■©'■^ “«(!)'<■ 


,j(0 '[(1 + »")% 4„^] + 2D^(^'( 1 + „^) + lDg^( 1 + 


3. /n = A:.^l,n=/ = l 


11 _ 3r, / 

^mlml o^lll 


.4(D„.2d4^^(0^ '(!+”■")+ 8022©“ 


12 21 [Sr, , 2 2 21 


“mlml - “mlml “ “ “^22 I ^1^12 ^^66 




»“l»l = |D„©"[(l+mV + 4»."] + 2D, 


© (1+m^) + lDg^(l+m^) 


©tal = %1„1 = 2(D„ + D«,)0(f)(l+».1 




4. m = k7il,n=l7il 


Q [(1+m^) +4m^] + 2(Di2 +20^6)0 (f) (1 + m^)( 1 + n^) + ^ 22 ( 1 ) [(1+n^) + 4n^" 
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= a 


= -{D„(f)"[(l+».V + 4»."] + (D,2 + 2D,,)Q(D"(l+m")(l+„")| 

“I™ = %L„ = -{D22(f)'[(l+"V + '>n'] + (Ol2 + 2D«)0'(D(l+m')(l+»")| (B-6) 

= Ol,Q'[(l+mV + 4m"] + D^(2)'(l+m")(l+n') + De,(l+m") 

= (D,2 + 0^)©(f)(l+'«')(l+«') 

4L« = £»22(f)^[(l+'^V + 4n^]+2D66QV+w^)(l+«^)+£»2y(l+«^) 

5. m = k = l,n-l = 2 

‘■ill, =-[‘>O„Q'42(D,,.2D«)0 ^(|)\i.O^|d45)\i.o"] 

»Sll = »?il, = + 

“Zu = “Jil, = |oJf)’(l+;)* + (D,2 + 2D,,)Q'(D(l+0" (B-7) 

HnU = -[‘'Dli©' + D66(f)"<l+')' + DgJ 

»f,l, = of,!, = -(Ol2 + D66)Q(f)<l+0' 

“Zu = -[|D22(f)'(l 4')'4^>66(f)'(l V 

6. m = k^l,n-l = 2 

‘■iL; = -^|o„(f)'‘[(l+».')" + 4m"] + 2(D,j + 2D,,)Q'(D'(l+m")(l+0^D22(f)‘'(l + 0‘'| 

»»»„; = »»L(= ^{D,i(f)"[(l+».V + 4».-] + (Di2 + 2D,,)0(0"(l+m")(l+/)"| 

“iL/ = %L/ = ^{D22(0'(l+O‘' + (D,2 + 2D«)0'(0(l+m')(l+O'| (B-8) 


24 



a 


a 


a 


22 

23 

33 
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a 


a 


a 


11 

mlkl 

12 

mlkl 

13 

mlkl 

22 

mlkl 

23 

mlkl 

33 

mlkl 

8 . 


a 


a 


a 


a 


a 


a 


11 

12 

13 

22 

23 

33 


m-k = 2,n=l=l 

= %‘«i = |o„Q u+«‘'+(o,2+2o^)Q(f)©+«" 

= %1H = '‘Dj2(f)’ + (O,2 + 2D„)0'(D(l+t)" (B-9) 

= -[|°ii(f)©V‘ 

= %ltl =-(0,2 + 0«)Q(f)(l+«' 

= -h22(f)'-^M!)©**>'-^°eJ 

m- k = 2, n = 

= ©L = ^[Ol,(f)©+«‘' + (0,2 + 2Ds6)Q(f)'(l+t/(l+»')] 

= = ^[022(f)’[(l+"V + 4n'] + (0,2 + 2D66)©'(f)(l+«'(l+"')] (B-10) 

“ ■ l[°"©°«(0© + ‘)"(1 + "") + Ofi.l 1 + ‘)'] 

= %© = -|<0l2 + 0<*)©©<l+«'(l+"') 

' ■ l{°22(f) '[(1 + »")© 4n'] + D«,(^ '(1 + t)"( 1 + n") + Dg^C 1 + »') | 
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9. m-k = 2,11-1 = 2 


“Ikl = +«‘' + 2(Dij + 2D,,)Q'(D^1 +tAl +/)“] 

»L, = “lu = -j[o„Q’(l+*)‘' + (D,2 + 2O^)©(f)'(l+«'(l+0'] 

»™« = »™H = -j[D22(f)V+/)‘'+(Di2 + 2D^)Q"(5)(l+»:)"(l + 0"] (B-11) 

“Zu = i["ll(f)'< > + *>" + °66(f)'< > + ' + ')' + >+*>'] 

»™« = »™H = i(Di2 + D66)©(f)(l+«"(! + ')" 

»L, = |[o22(f)'(l+0'‘ + D66(f)'(l+t)'(l+/)^D2,(l + 0'] 
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APPENDIX C 

TABLE OF SPECIAL INTEGRALS 


Different types of integrals are needed to earry out the integration of strain energy expression assoeiated with 
the dome-shaped temperature profile heating. 

Special Integrals Needed for 4S Case: 

mjix kjix, a , 

cos -eos- dx = - ■, m = k 

Q a a 2 

f cos^^cos^-^dy = - ; n = l 

Jo b b 2 

r%inBfI^sin^sin^dx = -^ ---- ; m±k = cvcn 

Jo a a a ^[(^fn + k) -l][{m-k) -1] 


r\in^sm^-^sm^dy = - 
Jo b b b 


^^ _ 


; n±l = even 


Special Integrals Needed for 4C Case: 


/„ 


^ ■ Jix 2 jix ■ mjix ■ kjix , amk f 
sm — cos —sm-sm- dx = -1- 


^ \[{m + kf-l][{m-kf-l] [{m + kf - 9][{m - kf - 9] 


fo 


“ ■ 2 ;rx Jtx ■ niJtx knx , am\ (w -k )-\ 

sm —cos —sm-cos- dx = —1 - - - 


{m -k)-9 


+ kf-l ][(m -kf-l} [(m -i- kf - 9][(m - kf - 9] 


fo 


“ ■ 2 JTX Jtx mjtx ■ kjtx , ak\ {m -k ) + \ 

sm —cos —cos-sm- dx = - —1 -:- 


{m -k )+9 


M[(m + kf-\ ][(m -kf-\] [{m + kf - 9][(m - kf - 9] J 

m±k = even 
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fo 


“ ■ 3 Jtx niJtx kjix, 3a I {m +k)-I 

sin —cos -cos- dx = - —■!-- 

a a a 


M[(m + kf-\ ][(m -kf-\] [{m + kf -9][{m -kf-9] 


■ 3 Jtx ■ mjix ■ kjtx , 3amk\ 1 1 

I sin —sin-sin- dx = -J--------- 

Jo a a a Jt [[(rn + k) -\][{m - k) -[{m + k) - 9][{m - k) - 9] 


Ain^cos^ ^sin^sin^dy = ^ 

Jo b b b b ^ + 


-!][(«-/) -1] [(« + 0 -9][(n-/r-9]J 

n±l = even 


(P ■ 1 Jty Jiy ■ mjiy kjiy , bn\ {n -I )-I {n -I )-9 

f sin -^cos-^sin—^cos-^dy = —J--------- 

Jo ^ ^ Z’ b + [{n + if - 9\{n - if - 9} 


(P ■ 1 Jty Jty njty ■ IJty , bl \ («-/) + ! {n -I ) + 9 

f sin ^cos ^cos^sin^dy = - —J-^---^- - - 

Jo b b b b [(„ + ;)2 _ 9 ]|-(„ _ ;)2 _ 9 ] 


/„ 


P ■ 3 Jty tiJty IJty, 3b { {n + 1 )- I 
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APPENDIX D 
BUCKLING EQUATIONS 

Case 1: 4S condition 

m±n = even (symmetric buckling) 


^kl~* ^11 "^13 ^22 ^31 "^15 ^24 ^33 ^42 ^51 "^35 ^44 ^53 



(D-1) 



Case 1: 4S condition 


m±n = odd (antisymmetric buckling) 



(D-2) 





Case 2: 4C condition 


m±n = even (symmetricbuckling) 

^11 ^13 "^22 ^31 "^15 "^24 ^33 "^42 ^51 "^35 "^44 ^53 



(D-3) 



Case 2: 4C condition 



(D-4) 
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FIGURES 




Figure 2. Shifted dome-shaped temperature profile heating (T ^ 0). 
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(a) 4S fixed. 



Figure 8. Two types of boundary eonditions for simply-supported edges. 
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(b) 4C free. 


Figure 9. Two types of boundary eonditions for elamped edges. 
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(a) Uniform temperature profile heating. 



(b) Dome-shaped temperature profile heating. 



Figure 10. Buekled shapes of square panel {a/b = 1) under different temperature profile heating; 4S fixed edge 
supporteonditions; = 0. 
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(a) Uniform temperature profile heating. 



(b) Dome-shaped temperature profile heating. 



(c) Roof-shaped temperature profile heating. 


Figure 11. Buckled shapes of square panel {a/b = 1) under different temperature profile heating; 4C fixed edge 
support conditions; = 0. 
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Figure 12. Plots of buckling temperatures functions of boundary heat sink temperature TJT^ for different 

temperature profile heating; fixed cases. 



Figure 13. Plots of buckling temperatures (7’^)p^as functions of boundary heat sink temperature TJT^ for different 
temperature profile heating; free cases. 
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Figure 14. Plots of buckling temperature magnification factors of the first and the second kinds { 17 , as 
functions of boundary heat sink temperature TJT^. 
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Figure. 15. In-plane deformed shape of square panel under different temperature profile heating; four edges fixed; 
u{x,y) 0, v{x,y) 0 at panel interior; half panel plot. 


43 













Figure 16. Distribution of x-displacement u{x, y) in square panel {a/b = 1) under dome-shaped temperature 
profile heating; T^ = \ °F; 4S fixed; u{x,y) 0, v(x,>') 0 at panel interior. 
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Figure 17. Distribution of ^-displacement v(x,y) in square panel (a/b - 1) under dome-shaped temperature 
profile heating; = I °F; 4S fixed; u{x,y) 0, v(x,>') 0 at panel interior. 
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Figure 19. Distribution of y 
heating; T^ = 1 °F; 4S fixt 
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Figure 20. Distribution of o^{x,y) in square panel (a/b = 1) under dome-shaped temperature profile heating; 
= 1 °F; 4S fixed; u{x,y) 0, v{x,y) 0 at panel interior. 
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Figure 24. Distribution of Oy{x,y) in square panel (a/b = 1) under roof-shaped temperature profile heating; 
Tg = 1 °F; 4S fixed; u{x,y) 0, v{x,y) ^ 0 at panel interior. 



Figure 25. Distribution of T (x,y) in square panel (a/b = 1) under roof-shaped temperature profile heating 
Tg = 1 °F; 4S fixed; u(x,y) ^ 0, v(x,y) ^ 0 at panel interior. 



Figure 26. Distribution of o^{x,y) in square panel (a/b = 1) under dome-shaped temperature profile heating; 
= 1 °F; 4S fixed; u{x,y) = v{x,y) = 0 everywhere. 



Figure 27. Distribution of Oy{x,y) in square panel (a/b - 1) under dome-shaped temperature profile heating; 
= 1 °F; 4S fixed; u{x,y) = v{x,y) = 0 everywhere. 
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Figure 28. Distribution of o^{x,y) in square panel (a/b = 1) under roof-shaped temperature profile heating; 
= 1 °F; 4S fixed; u{x,y) = v{x,y) = 0 everywhere. 



Figure 29. Distribution of o (x,y) in square panel (a/b = 1) under roof-shaped temperature profile heating; 
T = 1 °F; 4S fixed; u{x,y) = v{x,y) = 0 everywhere. 
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